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Abstract 

The sensitivity to initial conditions and relaxation dynamics of two-dimensional 
maps are analyzed at the edge of chaos, along the lines of nonextensive statistical 
mechanics. We verify the dual nature of the entropic index for the Henon map, one 
(Qsen < 1) related to its sensitivity to initial conditions properties, and the other, 
graining-dependent (q re i(W) > 1), related to its relaxation dynamics towards its 
stationary state attractor. We also corroborate a scaling law between these two 
indexes, previously found for z-logistic maps. Finally we perform a preliminary 
analysis of a linearized version of the Henon map (the smoothed Lozi map). We 
find that the sensitivity properties of all these z- logistic, Henon and Lozi maps are 
the same, q sen = 0.2445 . . . 
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1 Introduction 



The rebirth of the focus on nonlinear maps has opened unsuspected new- 
perspectives in physics. Among them, the logistic map [1,2] has achieved a 
paradigmatic rank. But not only dissipative maps have been playing impor- 
tant roles; also conservative and symplectic ones, like the standard map, have 
their places. These simple nonlinear models exhibit essential features of or- 
dered and chaotic systems, and also of the transition between them. When we 
deal exactly with the critical points, these models present many typical be- 
haviors of complex systems, for instance, slow dynamics, break of ergodicity, 
power-laws, self-similarity, among others. 

One of the important contributions that comes from the study of these maps 
aims to understand the fundamental connections between dynamics and ther- 
modynamics. In fact, the relation between these two main branches of physics 
is not yet fully understood and there are many points that deserve further 
investigations [3]. Once these systems have low- dimensionality, it is possible 
to computationally investigate them in its full phase space, something that is 
very hard (not to say impossible) to implement for high-dimensional systems. 
There are also important analytical results regarding the logistic maps at the 
edge of chaos [4,5]. 

It is very impressive that features exhibited by high-dimensional Hamilto- 
nian systems [6] may be found in some of these maps (conservative [7] or 
even dissipative [8]), e.g., the met a- stationary non-Boltzmannian state with 
a crossover to a stationary Boltzmannian state (crossover that is shifted in 
time as the system grows), indicating that the scenario conjectured in [9] is 
somehow ubiquitous. 

There are basically two computational paths for investigating such dissipa- 
tive systems, along nonextensive lines: sensitivity-based and relaxation-based 
approaches. By nonextensive we mean the use of the nonextensive entropy 
(written in its dimensionless form) 



and the concepts it implies. The first step of both approaches is to define 
a partition of the phase space into W equally spaced cells (the graining). 
Then the graining is made finer (by increasing W). In the sensitivity-based 
approach, we choose one of the W cells of the phase space, put within it 
N points (typically iV ^> W), and then follow the dynamical evolution of 
these points. Its initial entropy is S q (t = 0) = 0, Vg, (complete knowledge) 
so it increases with time. The chosen initial cell is taken as that (or those) 
that rapidly spread the points over the phase space. It is possible to estimate 
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the value of the entropic index, here called q sen , by at least three different 
methods, all of them giving the same results. To see a description of these 
methods, when they were first employed, see [10,11,12]. 

In the second (relaxation-based) approach [13], we initially take iV 3> W 
points uniformly distributed in the whole phase space. Its initial entropy 
S q (t = 0) = (W 1 ~ q — 1)/(1 — q) is a maximum, Vg (complete lack of in- 
formation). At the edge of chaos and as t — > oo, the distribution of points do 
not occupy the whole phase space, but a fractal subspace of it. So the entropy 
decreases (relaxes) with time. The shrinking of the occupied volume permits 
to estimate the proper entropy index, here called q re \. 

These two approaches to the same problem generate different entropic indexes: 
the sensitivity-based approach generates q sen < 1, while the relaxation-based 
approach generates q re \ > 1. In Ref. [14] it was found for the logistic map that 
the entropic index q rd depends on the graining W, and that there exists the 
following scaling-law relation: 

q rel (W ^ oo) - q rel (W) oc W~^. (2) 

In the present work we further investigate this scaling relation for the Henon 
map (preliminary results were given in [15]). The existence of a scaling-law be- 
tween sensitivity-based and relaxation (graining dependent)-based approaches 
for different systems, including high-dimensional ones, was suggested in [16]. 
We also perform a preliminary study of a linearized version of the Henon map, 
according to the sensitivity-based approach. 

Before going on, let us make a few comments about nomenclature. In early 
papers on nonextensive statistical mechanics, there were no awareness of the 
duality of the q indexes, so they were simply called q, or, sometimes, q*. The 
existence of the duality was shown in [13] and then commented in [17]. In 
that paper it was adopted the expression mixing in the sense of sensitivity to 
initial conditions, and equilibrium in the sense of relaxation to the stationary 
state attractor. So it was then employed the nomenclature q m i X and q eq . In 
some later papers on the subject, this nomenclature has gradually changed. 
In the present work we adopt q sen in the place of q m i X , and q Te \ in the place 
of q eq , due to the generality of the concepts. Such adaptation of concepts and 
nomenclature is natural in a still evolving branch of science. 



2 Sensitivity and relaxation for two-dimensional maps 

We analyze two examples of two-dimensional maps. The first one is the Henon 
map [18]: 
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xt+i = 1 - ax 2 t + y t , 
yt+i = bx t . 



(3) 



The parameter b recovers the logistic map when it is set to zero. If b = 1, 
the Henon map is conservative. We are particularly interested at small values 
of b (dissipative systems) at the edge of chaos (critical points at which the 
standard Lyapunov exponent changes sign). The Table 1 presents the critical 
values of the parameter a for different values of b, as well as the limits of the 
phase space in each case. 

Table 1 

Critical parameter a = a c for the Henon map, at the edge of chaos. x m i n , x max , y m i n 
and Umax represent the limits of the phase space along x and y axis, respectively. 



b 


a c 


%min 


%max 


Vmin 


Vmax 





1.40115518. .. 


-1 


1 








0.001 


1.39966671. . . 


-1.05 


1.05 


-0.0015 


0.0015 


0.01 


1.38637288. . . 


-1.05 


1.05 


-0.015 


0.015 


0.1 


1.26359565. . . 


-1.1 


1.1 


-0.10 


0.12 


0.2 


1.14904648. .. 


-1.179 


1.179 


-0.236 


0.236 



We also analyze the approximation of the Lozi map made at Ref. [19]. The 
Lozi map is a piecewise linearization of the Henon map. In [19] the authors 
have made a smoothness in the original Lozi formulation by inclusion of a 
third parameter e: 



(4) 



x t +i = 1 - aF e (x) + y t , 
y t +i = bx t , 

and F € (x) is defined as 

\x\ if \x\ > e, 

~ (5) 

^ g e (x) = (x 2 /2e) + (e/2) if \x\ < e, 
with < e < 1. The limit lim e ^ + F € (x) = \x\ recovers the original Lozi map. 

In order to establish a connection between q sen and q re i, we follow the proce- 
dure adopted in Ref. [14], that we briefly describe in the following: we begin 
with iV points inside a properly chosen single cell (similar to the sensitivity- 
based approach). The initial cell is chosen as that one which gives the maxi- 
mum spread of the points — a maximum overshoot in the curve S qsen vs. time 
(see Fig. 1(a), illustrated with the Henon map). The entropy is evaluated with 
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the value q sen , as determined by the sensitivity-based methods. Fig. 1(a) may 
be divided into two parts: the "uphill side" is governed by the sensitivity to 
the initial conditions; the "downhill side" is governed by the relaxation to the 
final stationary value of S qaen . These "two sides" are a visual expression of 
the non-commutability of the limits lim^oo linm/^oo and limyK^oo liim,-^. 
The former expresses that complete knowledge never meets a stationary state. 
On the other hand, the later order of limits expresses that a stationary state 
is only possible with lack of information. This non-commutability was con- 
jectured in Ref. [9] for many-body Hamiltonian systems, if we substitute the 
present variable W for the number of particles. 



We particularly focus on the "downhill side" and define the variable AS q = 



S. 



««»(*) - s q.e n (t -> °°). The time evolution of AS qsen gradually displays a 
power-law, as W (and proportionally N) increases (Fig. 1(b)). The slope (in a 
log-log plot) of the curve, in the range AS qsen follows a power-law, is identified 
with the power of a decreasing g-exponential (without cutoff), e q (—t) — [1 — 
(1 — q re i) t/T q ] 1 /( 1 ~ qrel ) (|slope| = l/(q re i(W) — 1)). It is somehow subjective to 
decide when the power-law begins and when it ends. So we average over many 
possible intervals (different beginnings and endings), taking into account those 
with higher correlation coefficient. 
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Fig. 1. (a) Time evolution of the g sen -entropy for the Henon map; (b) Time evolution 
of AS qsen . Different values of W, with b = 0.001, at the edge of chaos. 



The sensitivity-based approach to this map was applied previously [20] , when 
it was found q sen = 0.2445 . . ., the same value found for the logistic map 
[10,11,12]. When the relaxation-based method is applied to the Henon map 
[15, Fig. 1], once again it was found the same value of the logistic map for 
infinite fine graining, q re i{W — > oo) = 2.41 .... These results may include the 
Henon map into the family of the logistic map. We then examine the relation 
between these two q indexes: we approximately find the scaling law, Eq. (2). 
The scaling law is found when we analyze the projection of the map in its 
logistic subspace (x axis), so we use W x instead of W in Eq. (2). Fig. 2(a) shows 
our results for the Henon map. The errors associated with the obtained results 
are relatively large when compared with those of the logistic map. This may be 
due to the fact that two-dimensional maps demand higher number of cells W, 
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and consequently higher number of points N. Once the computer capabilities 
are, of course, limited, the graining we can achieve in two-dimensional maps 
are coarser than in one-dimensional ones. For this reason, the extrapolated 
curves for q re i(W — > oo) in Fig. 2(a) for b ^ should be viewed as indicative. 
We also conclude that the logistic limit b — > is slowly reached. We call 
attention that in Ref. [15] we were not able (at that moment) to see that it 
was necessary to analyze the logistic subspace, once we had only one value of 
b; in that paper we adopted W in the scaling law. Only now we amend that 
conclusion: Fig. 2(b) shows that the use of the full partition of the phase space 
W, in the two-dimensional Henon map is not fair, once b — > does not recover 
the logistic case. 
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Fig. 2. Scaling law for the Henon map (Eq. (2)), at the edge of chaos, (a) The 
abscissa is the partition of the logistic subspace, W x . (b) The abscissa is the full 
partition of the two-dimensional phase space, W, as appeared in Ref. [15]. The 
Figures make evident that the choice (b) is not fair, as explained in the text. 



Let us now address our results for the smoothed Lozi map, Eq. (4), at the 
edge of chaos (when e = 0.3, a c = 1.72879898... for b = 0.001, and a c = 
1.71996303... for b = 0.01). We have focused this map according to the 
sensitivity-based approach. The analysis of the sensitivity to initial conditions 
at the edge of chaos (see [10] for details of this method) shows a vanish- 
ing standard Lyapunov exponent. The positive g-Lyapunov parameter, \ qaen , 
expresses that the system follows a complex dynamics, whose upper bound 
of the sensitivity £(£) evolves through a slow dynamics, = [1 + (1 — 
(isen) K m t] 1 ^-*™), with £(t) = \im Ax{0) ^ Fig. 3(a) identifies the slope 
of the upper bound with the power of the previous g-exponential equation for 
£(£). A different method related to this sensitivity-based approach, namely the 
rate of increase of S q [12], is also presented (Fig. 3(b)). The proper value of q is 
that one which generates a linear increase of entropy (a finite rate of entropy 
production). Both methods give q sen ~ 0.2445, the same value of the logistic 
map. It is still missing the analysis of this smoothed Lozi map according to 
the relaxation-based method, as well as the correlation between its entropic 
indexes q sen and q rel . 
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Fig. 3. Smoothed Lozi map (Eq. (4)), with e = 0.3. (a) Time evolution of the 
sensitivity to initial conditions, (b) Time evolution of S q . 

3 Conclusions 



The duality among the g-indexes, one related to approaches based on sensitiv- 
ity to initial conditions, q sen , and other related to approaches based on relax- 
ation to stationary states, q re i, and the connection between them, being verified 
for the first time for the logistic map (Ref. [14]), has been approximately ver- 
ified in a two-dimensional (Henon) map, which belongs to the logistic family. 
The relation is verified in the logistic subspace of the two-dimensional map. 
The results straighten the validity of the scaling relation, Eq. (2). Possibly 
similar scaling laws exist for different phenomena, as suggested in Ref. [16], 
and further studies along these lines would be very welcome. The prelimi- 
nary analysis of the smoothed Lozi map, that introduces a linearization in the 
Henon map, indicates that, as far as the sensitivity to initial conditions are 
concerned, this map presents the same value of q sen as that of the logistic and 
Henon map. Further studies are necessary to evaluate the characteristic values 
q re i, as well as the connection between them. Recently it has been conjectured 
in [21] that the complete characterization of a system may demand not only 
two, but three entropic indexes. The third one would be related to the station- 
ary state itself, named q s tat- These indexes would certainly be related among 
themselves. Detection and estimation of q s tat, as well as the discovery of how 
these three indexes are connected, are essential steps towards the closure of 
the nonextensive statistical mechanics theory. 
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